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INTROmCTlON 

Hus investigation began with a desire to locate the Gulf Stream by detecting the 
geost rophically induced slope in the ocean surface across the current. The Gulf Stream 
was chosen because: (l) a 1 meter surface elevation change was expected across the Stream, 
(2) it was one of the most studieu currents, (3) it was in the <2 : 0S 3 calibration area, 

(4) a detailed geoid was available for the area, and ( 3 ) corroborating information defin 
mg the location of the Stream was available from the Coast Guard, NOAA-fCSS, and the 
Naval Oceanographic Office. If the Gulf Stream could be detected, then two other western 
boundary currents could also be studied; i.e., the Kuroshio Current south of .Japan ;unl the 
Agulhas Current east ol Africa. The analysis began by selecting a set of 20 northbound 
CTOS altimeter profiles, both global ;ind intensive modes, which crossed the Gulf Stream. 
Northbound passes were chosen because they would cross the current at a more orthogonal 
.ingle in the region between Florida and Cape Hatteras and thus provide a maximum slope 
when crossing the Stream. Both global and intensive mode data were u sod to evaluate the 
advantages of each mode. The early analysis showed that the intensive mode data, with its 
inherently lower noise characteristics required less filtering ;ind therefore shorter 
length filters could be used. Also, the 100 sample per second data was not used for three 
reasons: (1) it was felt that 10 samples per second data was dense enough to detect the 
elevation change, (2) much longer filters would be required to force the 100 samples per 
second data to provide results compatible with the 10 per second data, and (31 not much 
100 sample per second data had been acquired. 



\ tew di He rent altimeter measurements and der i ved parameters were investigated to 
determine it thev produced am .inomalv in the signature when sensing the (lulf Stream. 

Hie altimeter measurements studied were sea surface heights (SSI! IT!- ) , automatic gain 
control (AfiC) # and the altimeter waveform .ample .tnii hold gates. In addition, the derived 
parameter ^ea state (H 1/3) wa- also studied. In investigating the SSHITT capability, 
correlation of both the geo id .ind ocean bottom topography with measured altitude was 
evaluated. Preliminary analysis indicated that the SSIUT! measurement when referenced to 
a detailed gee i d showed the most promise in being able to detect the flulf Stre.un. The 
other parameters, A (IT, H 1/3 ind the waveform data did not show much potential at that 
time, so that total effort was directed to the S9UTI analysis. Incidentally, this tech 
nique had been tries! once bet ore (ref 1) using SKY LAB altimeter data referen ed to a 
15" x 15” geo id. Hie results were inconclusive. 

DATA ANALYSIS PP0CTI1JRE 

A satellite radar altimeter is an instniment which measures precisely the distance 
from a satellite platform to the ocean surface. The geometry of the measuring system is 
illustrated in figure 1. The altitude measured by the altimeter can be represented as: 

h a “ h s ' V Ah * ‘n (1) 

where is the altitude measured by the altimeter, h^ is the satellite height above a 

reference ellipsoid as estimated from satellite tracking data, is the geoid height, 

Ah is the height deviation from the geoid due to oce.in dynamic processes such as tides 

.ind currents, and • is the random error in the measurement. Initial analysis did not 
n 

consider Mde effects because of the short wavelength of the profiles being analyzed. 
However, later on in the analysis longer profiles were analyzed, thus requiring a consid- 
eration of tide effects. Typical ranges and uncertainties of the terms in equation (1) 
are shown in Table I . 

Tne altimeter was designed to make 100 individual measurements per second and to 
telemeter either them, or a 10 pulse per second average, to the telemetry stations. For 
the reasons mentioned above and to save processing time, only the 10 samples per second 
average is used. Ihese altitude measurements are then preprocessed and converted to 
sea surface heights (ref. 2) and are treated in this report is raw sea surface heights. 

The data processing procedure, using data from orbit 1710, is sumarized in figure 2. The 
raw sea surface heights shown in figure 2a are first edited to e: iminate anomalies due to 
internal noise and unknown data spikes. T l\e edit criterion is based on a predicted sea 
surface height calculated by fitting a straight line through the last eight seconds of 
data (80 data jx^int^ representing 60 km in physical distance). .Any point differing from 
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the predicted! height hv more than two meters ( approx imat el v three standard «k»viations of 
the noise level of ’’O an) is replaced hv the predicted height value, Ihe edited sea 
surface height data shown in figure 2h was then originally filtered *'ith a 41 f>oint, equal 
weight, midpoint filter. Ihis filter did not suppress enough noise end resulted in very 
noisy velocity estimates. After a brief filter study, .in HI point equal weight, midpoint 
filter was chosen since it reduced the noise without seriously compromising the sea sur 
face signature and resulted in mor? stable velocity computations. Hie eight second 
time constant was chosen so that tne noise level of the data could be maintained below the 
10 cm level. After the filtering process, the smooth sea surface height was referenced 
to the ‘tarsh-Giang S’ x 5* geoid (see figure . This geoid, however, was onlv defined 
between lb° N to 39° N .ind 2* ? 8° F to 300° F. Ye , the Gulf Stream could verv r easily 
meander north of the latitude limit. Therefore, mean anomaly gravity data was obtained 
and a new S’ x S’ geoid was computed to 41.5° N thus providing mnre geoid definition in 
the Gulf Stream meandering area. Subtracting the geoid from the smooth sea surface heights 
results in a residual which is nearly flat in the static open ocean, as shown in figure 2d. 
Next, to minimize the error between the geoid and the smooth sea surface height south or 
cast of the mean position of the Gulf Stream, a linear fit is made to the residuals over 
the section representing the open ocean. The straight line is then subtracted from all 
the residuals thereby removing any potential orbital bias or slope errors .uni producing 
an estimate of dynamic heights as shown in figure 2e. Finally, the dynamic heights are 
differentiated with a filter designed to match the smoothing fil er in order to compute 
the slope of the sea surface along the profile. The computet! slope is substituted in 
the following geostrophic equation to obtain velocity estimates as shown in figure 2f. 
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& Ah 

sin * aT 
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where is the surface velocity; g is the gravity acceleration at 980 cm/sec ; P is the 
angular speed of the earth at .0000729 rad/sec; * is the latitude of the location; Ah 
the total height anomaly ;ind aL the horizontal distance over which the height anomaly 
occurred. 


RESULTS 

'Hie first few months of this investigation were devoted to processing as many 
northbound passes as possible in the area of the Gulf Stream using the data processing 

procedure outlined above. Each region where the Gulf Stream flowed had its own problem 

areas. For example, in the region east of the Florida and Georgia coasts there was not 
enough altimeter data acquired over the open ocean to supply enough data for a reliable 

fit. Either the profile would intersect land in the Bahamas or the altimeter was not 



3 


turned on until after passing over the Isl.md* . Also , the accuracy of the Marsh ()i;ing 
geo id in th<s area was suspect. of the earlv analysis was performed in the area 

he tween South Carol ma .uni G.i|x* Hatteras. Iliis area was chosen because the suhsatellite 
trace crossed the G*ul f Stream in .ui almost orthogonal angle .uni there ton' would provide a 
nuximum slojx* signature in the profile. There were two shortcomings in this area, f irst, 
the current sterns to hug the break in the continental shelf. Thi^ characteristic made 
both the geo id .uni the Gulf St re. im produce .ui anonulv along t ho same general area, making 
earh serration of the two signals difficult. Also, the geo id near the Blake Spin was 
not defined well enough to produce the kind of results expected as can be seen in figure 
Note the double hump in the dv runic height, this occurs directlv over ui ocean bottom 
feature associated with a steep gravity gradient. Next , the region centered around i ape 
Hatteras provided early difficulties also because of the sharp geo id ch.inee occurring 
along the line of current flow. 1 n addition, the geoul accuracy was igain uncertain. 
Finally, the analysis moved out to the open ocean east of Cape Hatteras where it was felt 
at the bottom topographic influence on the Gulf Stream would be minimal uni that the 
geo 1 d , especially between J8S° I .uni 292° 1 would be relatively flat. The analysis in 
this region provitled the first conclusive proof that indeed the slope of the ocean surface 
was located in the s.ime vicinity win c the NO.AA satellite 1R imagery indicated .i surface 
temperature anomaly indicative of a thermal front. Analysis in this area progressed unt i 1 
it was found that profiles in the vicinitv of the Vw Ingl.unl Seamount Cham tended to 
produce jxxm results. Hie jxx>r distribution ;uui lack of measurements in the area of the 
se.imounts pnxluced large uncertainties in the wan anomalies m this general region. 
Therefore, profil * were terminated earlv or started late when approaching the seamount 
area. 


Ocean Truth Sources 

In order to be certain that the altimeter profile actually passed over the Gulf 
St re.im, s t irx' source of ocean tmth (Gulf Stream location) was neeiled. The first product 
used was that of the "Gulf Stream,” a monthly publication produced bv SOW . ihe in forma 
t ion in it is good, however, usually onlv two western boundaries are given; one for the 
beginning of the month and one for the end of the month. Hie publication is verv 
useful for following the path ot some spawned eddies which arc not followed bv other 
data products. The next data product used was the F.xper imental Gulf Stream Analysis 
produced weekly bv NOAA-NFSS. The product defines current width and is based primarily 
on satellite VHRR-IR imagery and therefore is limited bv cloud cover. The cloud cover 
especially in the area south of Gape Hatteras precluded some earlv analysis in this area. 

\ final product, the lxperimcnt.il Ocean Frontal Analysis produced hi -week lv has the 
advantage of supplementing the IR data with hydrographic data. This data product also 
was able to discern both the eastern and western boundaries of the Gulf Stream more readily 
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than .uiv of the other data product^ u*ed thus far. Therefore, roost analysis conjxirison 
were made using the fxper iroental <Ve.»n frontal Analysis to define Gulf Stream current 
positions. 


Single Profiles 

After the dvnamic heights wer«* confuted the profiles were analyzed to iletermine 
the location of the f»ulf Stre.im. The inferred eastern and western boundaries of the 
current were defined at the breaks of the sharp height changes of the dynamic height 
profiles. The altimeter determined boundaries were then compared with the boundary* de- 
fined by NOAA's Experimental Gulf Stream Analysis product . An exair^>le of a single pass 
analysis is shown in figure 4. Note the agreement of the altimeter determined western 
boundary and the NOVA -determined boundary. In order to be able to compare more altimeter 
passes with Gulf Stre.im truth data, and because the Naval Oceanographic office's anal vs is 
seemed more detailed and reliable, the Experimental Ocean frontal Analysis (EOFA) product 
became the independent baseline for comparing future altimeter inferred Gulf Stream 
boundary data. The mean position discrepancy m the determination of the western bound 
arv bv the altimeter and EOFA was calculated for bO profiles acquired during August, Sep- 
tember, and October 19^5. Hie bias of the al t imeter* determined western bouidary was 2 A km 
♦ 3b km west or north of the highest surface thermal gradient as defined bv the EOFA chart 
Hus result is encouraging since Hansen and Maul (ref 3) determined that the thermal sur 
face boundary was 14.5 ♦ 11.8 km north or west of the Gulf Stream axis defined to be 
located at the position of the 15° C isotherm at 200 m depth. This results in the 
altimeter western boundarv being approximately 3" km west of the core of the current, a 
reasonable result. In addition, .in estimate of the height difference and surface velocitv 
was made us^ng the s.une bO passes. The mean surface current velocitv was 107 ♦ 29 cm/sec 
and the mean height difference was 100 ♦ 15 cm. Both of these means are in very much 
agreement with both measurements of the surface velocity .anil with expected height devia- 
tions across the Stream. 

Initial results were further substantiated by con^aring the results of consecutive 
passes occurring over the same profile. Results of three profiles acquired over the s.aroo 
track over a three month period are shown in figure 5. Note the agreement of t ho signa- 
tures of both the dynamic topography and the derived surface velocity. Hie relative agree 
ment of these passes was verv encouraging since the decreasing surface velocity trend was 
in agreement with reported mean values (ref. 4). 

Another short stud;* was performed with data off the coast of Onslow Bay, North 
Carolina, acquired during the first week of Mav \9"S. Detailed ocean truth data was 
available from R. Perchal of the Naval Oceanographic office and it was decided to compare 
the two data sources. The results of this analysis are shown in figures b and Note 
the signature of the altimeter data compared to the special EOFA chart. It can be seen 


5 


1 \ ' 



from the chart that water is moving north at "C" ; however, it "B" a south component of the 
loop current or me angler h seen which then flow'' back in a north direction at "A". Looking 
at figure P it can be seen that points A, H, and C on the map correspond quite well with 
the changes in siojx* in figure Ihe total meander is shown to be ab<xjt ISn km wide with 
surface velocities on the order of ISO cm/sec. 

rhroughout all the analysis thus far, the only valuable "truth” data used to 
corroborate the altimeter results were the IOTA products based primarily on \1!RR IR 
satellite data and nindom ship reports anil oceanographic data. It was decided that in 
situ oceanographic data would be necessary to further justify the capability of altimetric 
techniques. Therefore, an experiment designed by N. 1 . Huang was lr^lemented during May 
;ind lune of lSTb to compare hydrographic data with remote sensing measurements from in 
strunentat ion such as satellite radar altimeters and aircraft laser prof 1 lometers and a 
radar scatterometer. A result ohtainevl over the four day period 19 May to 22 May 19?t> 
is -’.own in figure 8. The darkened contour is the boundary' of the Gulf Stream as deter- 
mined from the 1:0 FA charts. Two tongues of water are discerned; one east of South 
Carolina ;ind the other east of Cape Hatteras. Passing through each of these tongues .and 
the (iulf Stream is a GHOS altimr' r pass. A plot depicting the estimated surface \elocit\ 
across the subsatellite track is shown for two profiles. In addition, .another profile 
crosses the Stream in the southbound direction. A ship crossing the Gulf Stream from Cape 
Charles, Virginia, was acquiring hydrographic data from which surface velocity could be 
esti.nated. Hus velocity i also [ lotted in figure 8. Note the illative agreement of the 
sij. ace velocity estimates which intersect the current perpendicular to the flow. The 
agreement of t ne hydrographic .and altimetric profiles provides the most dramatic support 
to the altimeter data. In addition, the velocity estimates across the tongues of water 
being shed from or coalesced with the Gulf Stream are indicative of the movement of these 
waters, figure 8 was obtained from a detailed hydrographic analysis on this two-week 
experiment prepared by T. B. Curtin of North Carolina State University. 

Dynamic Topographic 'laps 

All results so far have shown only instantaneous estimates of ocean surface dynamics. 
Maps attesting to deline tc the posit ion of the Gulf Stream were prepared using onlv 
northbound data acquired during August, September, and October lO^S. It mist be pointed 
out that these mips are not monthly means, but rather a mosaic consisting of instantaneous 
measurements obtained at different times during the month similar to a scanning time 
exposure along the subsatellite trace. The direction, number and both time and spatial 
distribution of the satellite passes influence the validity of the maps to a large degree. 
The contour mips were based on dvnamic topographies contoured at the 2n cm level using 
data located east of Cape Hatteras and avoiding the New Ingland Seamount Chain. A 



typical result i*> shown in figure 9 computed for the month of September. rhe crowding of 
the contour lines running northeast from tape Matteras agree quite well witn the general 
direction ot flow of the Gulf Stream as shown in reference 5. Some shortcomings of the 
above method are that it did not use .uiv southbound data passes, it did not use am passe** 
intersecting land south of Cape Matteras and that the fin.il product was not truly a mc.in 
map. Hiese problems were overcow by enlarging the data base to include southbound tracks 
ajul tracks intersecting land south of (ape Matteras. The southbound data was used to 
constrain the northbound passes by using a mininum variance technique to minimize the 
differences of the profiles at the intersecting ixmit**. This tec’inique is documented in 
reference t> which describes the theory and opera** n of the conputer program, si MIT. 

Results obtained f rum employing the technique are aown in figure 10. In general, the maps 
were obtained b' processing the data acquired ov< one month through the SI.AIfl program. 

Hie adjusted value's obtained at the intersecting ,>oints were the point* used to bo contour 
ed. \oto the absence of valid contour lines southeast of orbit 309". .he absence of 
intersecting points in the southeast area caused this void to occur. 

In a further extension of processing the data acquired during the latter six months 
of 1 9"5 , it was decided to calculate a mean surface dvnamic topography over the entire 
six months of data. The mean surface was confuted in the following manner. First, all 
six months of dynamic height profiles were analyzed through the SI MU program. Next , 
in places where many profiles crossed one another, a cluster was defined (see figure 11). 

\ mean and standard deviation ot the values in the cluster were calculated using the 
adjusted values from the SKMIT program. Hie cluster means over the six month period were 
contoured os shown in figure 12. The pass coverage used to produce the map is in the 
bottom right of the figure. This map shows a well defined Gulf Stre;im initially flowing 
north from Florida ,uui then changing direction about 3J° \ ,uul then flowing in a north 
easterly direction. Vlso note the depression of 40 cm at location 34° N, T’ W. This 
low was caused bv im eddy maintaining a static position over the span of three consecutive 
months during the six-month data span. Also, the high elevation in the northeast comer 
of the map is indicative of the New Ingland Seamount Chain influencing .he results since 
the geoid is poorly defined in that area. In addition to the mean topography over the 
**:x month period, the standard deviation is shown in figure 13. This map essentially 
depicts the regions where high energy cells exist in the area ot computation. For example, 
in the region north of 35° N, large anomalies exist where the Gulf Stream meanders quite 
freelv. Smaller isolated depressions exist in the central portion of the map. These 
depressions tend to occur where the cold eddies, having spawned from the current, are 
migrating in a southwesterly direction. Finally, near the coast of Long Ray, North 
Carolina, another dynamic region is apparent. Hus coincides with the meandering at the 
coast which has been mentioned (ref. " and 81. Next, the map for the monthlv means was 
calculated using onlv the values for that particular month which were available within 



the cluster. A typical plot is shown for November in fi>!iire 14. \gain note the definition 
of the Gulf Stream ;tnd the mr.imter north of 35° \ between 70° W ;ind 65° W. 1 1 rut 1 1 v # a 
Sv*a surface dynamic height anomaly is calculated by differencing the mean values for each 
month with the overall six month mean. A representative plot is shown in figure 15 for 
the month of Severer. The neander above 35 c N is easily seen in this type presentation 
a* are the eddies just soutn at 55° * 

A dramatic application of using the sea surf ** dynamic height anomaly maps can 
he seen in figure lb. The tracking of a Gulf Stream ring using the depress nxis found in 
the sea surface with altimetrv is shown compared to the inferred position determined 
from satellite IR imagery. The observed six-month sequence shows that the agreement is 
verv goc-d in months when both techniques detect a ring. 1W> interesting points are 
apparent in this figure; first, the eddy seems to be stationary during August through 
October (both sensors agree with this conclusion) and second, during November t’.c 
satellite IR does not detect the ring probably because of cloud cover or a weak surface 
thermal gradient. However the altimeter, which is an all-weather system, continues to 
track the ring. 

All maps generated to date have been contoured manually. A simple technique to 
contour November data using a computer method is shown in figure 17. TV technique divides 
the scanned area into 5' x 5* cells. All adjusted observations occurring within a cell 
are averaged. Any cells with ^ observations are filled by using a linear interpolat ion 
along a coastant band of latitude. Some strange signatures exist using this method such 
as the anomalies at 30° N and 291° E and at 35° N and 29b° E. These few anomalies are 
due to the type of interpolat ion being employed. A new three-point planar method is 
being progranrved at present and it is hoped that it will resolve some of the strange 
signatures . 

Even though most work done during the execution of this invest igat ion was done with 
the Gulf Stream data, another area has been studied. As defined earlier, in order to 
employ this current detection technique, one must have access to a well-defined, high 
resolution geoid. Such a geoid was not available either off the east coast of Africa 
to study the Agulhas and Somali Current, or off the south coast of Japan to study the 
Kuroshio Current. Gravity data was available in the Japan area and these point anomalies 
were averaged to produce 5* x 5' mean anomalies. The mean anomaly data was then used by 
Detailed Geoid Computation Program (DGC) to compute a S' x 5' resolution geoid for the 
Japan area. Early results in this area are very poor primarily because: (1) The area of 
definition is centered on three trenches south of Japan and the steep gravity gradients 
are difficult to map with available spatial resolution, and (2) in some areas there 
existed almost no data and therefore bathvmetric data was used to interpolate and infer 
gravity data so that a geoid could he computed. It is hoped that more dense gravity data 



Will become available in the Japan area so that a better >*c % o 1 J can be compute: an*! there 
fore work can continue on the kuroshio Current. 

OWXUSIONS 


It has been demonstrated that the larger geos trophic ocean currents, their meander* 
.uni spawned eddies can be detected and measured utilizing the techniques employed in 
this study. Manv applications are apparent in the areas of weather forecasting, climate 
changes, submarine warfare, ocean fishing, and shipping industries. Radar altimctrv is 
onlv a surface remote sensor; however, it con provide the field of dynamic oceanography 
with sea surface parameters which up to the present were not directlv measurable. Since 
no measurements c.ui be made below the surface with altimetry, it should complement more 
traditional hvdrographic methods. 
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Figure 3. Dynamic heights profile and bottom topography for repeated orbits 1625, 2151, and 26 
depicting gcoidal error near the Blake Spur. 




Single pass profile for orbit 53(1 indicating the Gulf Stream'- eastern 
'Mm and western MV' boundaries as estimated from thi- studv ind the 
western boundary i- estimated by Vi\A-M S\. 
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Figure 5. Comparison of Dynamic Topography and Derived Velocity for three 
repeated altimeter profiles acquired over the same ground track 
during three successive months. 



Figure b. Composite drawing of IR data from NOAA-4 satellite and 
sea surface temperatures (SST) data from passing ships 
with lie If Stream boundaries delineated. Suhsatellite 
track for orbit 374 is superimposed showing the 
location of points \BfiC (see figure 7) . 
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Figure 8 . Composite map showing IK derived Uilf Stream boundaries with 
ship .uid 1*05-3 tracks superimposed. (Hip courtesy of 
T. I'urtin'). 
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Figure 10. 


Sea surface heights for Novenbei 19^5 from anal>sis adjusting 
data for crossings from only those orbits shown by minerals. 







DATA CLUSTER 



Figure 11. lixample of a data cluster used in analysis 
for ft months data. 
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Standard deviation of sea surface topographv estimated at each luster 
of intersections as shown ii inset, July -December 19"'. 


Figure 13. 
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Figure 14. Mean sea surface topographv using elevations retained at each crossing 
shown in inset, November 19“*5. 




Figure IS. Differences between overall and Novemhet ID'S "xan surface*. 






MEAN SEA SURFACE 



Figure l* 7 . Example of conputer plot of November 19”5 data processed 
using ill points f ran those orhits shown in figure 10. 




